Adiabatic multicritical quantum quenches: Continuously varying exponents 
depending on the direction of quenching 
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We study adiabatic quantum quenches across a quantum multicritical point (MCP) using a 
quenching scheme that enables the system to hit the MCP along different paths. We show that 
the power-law scaling of the defect density with the rate of driving depends non-trivially on the 
path, i.e., the exponent varies continuously with the parameter a that defines the path, up to a 
critical value ct = ac\ on the other hand for a > Oc, the scaling exponent saturates to a constant 
value. We show that dynamically generated and path (a) -dependent effective critical exponents as- 
sociated with the quasicritical points lying close to the MCP (on the ferromagnetic side), where 
the energy-gap is minimum, lead to this continuously varying exponent. The scaling relations are 
established using the integrable transverse XY spin chain and generalized to a MCP associated 
with a d-dimensional quantum many-body systems (not reducible to two-level systems) using adi- 
abatic perturbation theory. We also calculate the effective path- dependent dimensional shift do[a) 
(or the shift in center of the impulse region) that appears in the scaling relation for special paths 
lying entirely in the paramagnetic phase. Numerically obtained results are in good agreement with 
analytical predictions. 
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Following the Kibble-Zurek (KZ)[l[i| prediction of the 
scaling of the density of defect in the final state of a 
quantum many-body system following a slow quench0,|3| 
across a quantum critical point (QCP)[^ @, there has 
been an upsurge in theoretical studies on quantum 
quenching across critical points [7l-[l^(for a review see 
The KZ argument predicts that the scaling of the 
defect density (n) in the final state is universal and is 
given by n ^ \lf^d/{vz+i) -^yjigj-g j- jg ^j^g inverse rate 
of driving across a QCP with the correlation length and 
dynamical exponents v and z, respectively, and d is the 
spatial dimension. The possibility of the experimental 
verification of Kibble-Zurek scaling (KZS) in a spin-1 
Bose condensate [l^, in ions trapped in optical lattices 
pTl [isj , and also in ultracold fermionic atoms in optical 
lattices [l^, [l^l has paved the way for the above men- 
tioned theoretical studies. 

Although, the KZS for quenching through a quantum 
critical point is well-understood; the scaling of the defect 
density following an adiabatic quantum quench across a 
quantum multicritical point (MCP) is relatively less stud- 
ied. A non- KZS behavior (n - l/r^/^) of the density of 
defects (wrongly oriented spins) for quenching across the 
MCP of the spin- 1/2 transverse XY chain was reported 
for the first time in reference Q which was later explained 
in reference [2l| introducing an effective dynamical expo- 
nent Z2{~ 3) for Jordan- Wigner solvable spin chains [22| 
reducible to a collection of decoupled two-level systems in 
the Fourier space and applying Landau-Zener (LZ) tran- 
sition formula [2^. This argument was extended to the 
non-linear quenching of a general Hamiltonian in refer- 



ence |24l |. In a recent communication, Deng et aZ[25[, 
attributed this anomalous scaling behavior to the exis- 
tence of quasicritical points close to a MCP where the 
energy gap is minimum and proposed a generic scaling for 
the multicritical quantum quenches in terms of the effec- 
tive critical exponents associated with these quasicritical 
points. Studies on inhomogeneous quantum transitions 
across a MCP has also been reported recently [2^. But 
how does the scaling exponent depend on the direction of 
approaching the MCP7 In this work, we address this par- 
ticular question proposing a different quenching scheme 
that enables the system to cross the multicritical point 
along different directions and derive the corresponding 
KZS which reduces to all the previous results in appro- 
priate limits. 

At the outset, we summarize our main result using 
the example of a one-dimensional spin- 1/2 transverse XY 
Hamiltonian [23] given by 



(1) 
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where tr's are the Pauli spin matrices, Jx, Jy are the in- 
teractions along X and y, respectively, and h denotes the 
strength of the transverse field and henceforth, we set 
Jx-'r Jy=^- The parameter j = Jx — Jy is the measure of 
the anisotropy of the interactions and + Jy = 7 = 1 
refers to the transverse Ising limit [l^. The phase di- 
agram of the Hamiltonian is well-known and shown in 
Fig. 1; we study the quenching across the MCP (A) at 
/i = 1 , 7 = with the correlation length exponents and 
the dynamical exponents = 1/2 and zi = 2, respec- 
tively. 

We propose a generic quenching scheme where the pa- 
rameters h and 7 are both varied simultaneously follow- 



2 



ing the relation 

= 1 + |7(i)rsgn(t);a > and -f{t)^-t/i 



(2) 



with time t going from — oo to oo; the system reaches 
the MCP A at t = 0. The quenching path hits the MCP 
verticahy (i.e., \d^ / dh\^^Q — ^ oo) for a > 1 and hori- 
zontaUy for a < 1. The parameter a introduces a non- 
Unear temporal variation of /i, thereby determining the 
path of quenching across the MCP as shown in Fig. 1; 
for a = 1 (path I), the MCP is approached linearly and 
hence the path is always equidistant from the Ising and 
Anisotropic critical lines. Similarly, if the parameters 
h and 7 are changed simultaneously by the same non- 
linear rate as ^{t) = h{t) — 1 = |t/T|" sgn(t) as studied 
in [2^, the quenching path once again happens to be 
equidistant from these two critical lines. The same is not 
true for the quenching scheme (2) if a 7^ 1; paths with 
a — and a — > 00 correspond to the passage across 
the anisotropic and Ising critical lines, respectively. For 
a = 1, the scheme (2) is the same as that of reference 

Our studies reveal the existence of a critical value 
ac{= 2 for model (1)) for the quenching scheme (2). For 
1 < a < 2 , the defect density (712) scales as n2 '-^ r~"/^, 
i.e., the exponent varies continuously with the parameter 
a and saturates to the scaling 1/3 for a > 2. We argue 
that the quenching path crosses the quasicritical points 
for a < 2 but not for a > 2 when we recover the scaling 
relation valid for quenching the system across the MCP 
by linearly varying the parameter 7 along the Ising crit- 
ical line [Tjl that amounts to setting a — 00 in Eq. (2). 
We do also derive similar scaling relation for the special 
path III (Fig. 1) and calculate the path-dependent dimen- 
sional shift appearing in the scaling. Finally, we extend 
our studies to a d-dimensional generic Hamiltonian using 
the adiabatic perturbation theory 01 and establish the 



scaling relation, 712 ~ r 



-adu2 /[a{z2i^2 + 'i-)+ziV2['i--oi)] -^^fj^gj-g 



Zi^vi are the exponents associated with the MCP while 
2:2(0;) and V2{a) = viZi/ Z2{a) are path dependent expo- 
nents associated with the quasicritical points. To the best 
of our knowledge this path-dependent 22(0:) leading to a 
path-dependent continuously varying exponent of KZS 
were never reported before. It needs to be emphasized 
here that these path-dependent exponents are generally 
not useful in desrcribing the equilibrium quantum criti- 
cal behaviour though they do dictate the non-equilibrium 
dynamics of models as discussed above. 

It is well known that using the Jordan- Wigner trans- 
formation [2^, the Hamiltonian (1) can be reduced to 
direct product of decoupled 2x2 Hamiltonians for each 
momentum k given by 



h + cos k i"f sin k 
—i"fsmk — (/i + cosfc) 



and that the excitation gap is 



Afc 



'ih- 



cos fc)2 + 7^ ! 



(3) 
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FIG. 1: (color online) The phase diagram of a one-dimensional 
XY model in a transverse field. The vertical bold lines 
given by ft — ±1 denote Ising transitions from the ferro- 
magnetic phase to the paramagnetic phase. The horizontal 
bold line stands for anisotropic phase transition from a fer- 
romagnetic phase with ordering in the i-direction to a ferro- 
magnetic phase with ordering in the y direction. The mul- 
ticritical points are shown by points yl (ft = 1,7 = 0) and 
B{h — —1,7 — 0). We show different quenching paths corre- 
sponding to different values of a; path I (path II) is for a = 1 
{a — 2). Path III corresponds to a quenching scheme in which 
the system is always in the paramagnetic phase and touches 
the MCP at t = 0. We show that in the shaded region there 
is a continuously varying effective dynamical exponent and 
the exponent for defect density while for a > 2, n scales as 



A 

0.0006 
0.0004 



0.0002 - 





Path III / 






a = 0.9 /' 








a=l 















-0.0002 



0.0002 0.0004 0.0006 
Y 



FIG. 2: The variation of energy gap as a function of 7 for 
k = 0.01 close to MCP and also the quasicritical point for 
different q. 



In the vicinity of the quantum MCP A, with TT — k = k 
0, Eq. (3) can be simplified to 



iJfe = (/i-l + fc2)6-^+7A;CT^ 



(5) 



For the sake of comparison, let us briefly recall the 
results of quenching across the MCP A with j{t) = h{t) — 
1 = |i/r|"sgn(t) as studied in [2^ using the linearization 
process. The Hamiltonian (3) gets modified to 



|-rsgn(t) + fc2)a^ + |^psgn(t)M^ 
r r 
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(|-rsgn(t) + fc")<T"+fc^<T 

T 



(6) 



In deriving (6) an appropriate unitary transformation 
is used so that the time-dependence of the off-diagonal 
terms is removed and the LZ transition formula is directly 
applicable for a = 1 [1^ . We observe that the energy gap 
is minimum at a time so that |to/''"|"sgn(t) + fc^ = 0; 
this defines a quasicritical point at to on the ferromag- 
netic side of the MCP (i.e., t < 0). The minimum gap 
scales as so that we have an effective dynamical expo- 
nent Z2 = <i even though the dynamical exponent at the 
MCP (i.e., t = 0) is still given by zi = 2; the equivalent 
correlation length exponent 1^2 — (y\'^\)lz2 = 1/3 . We 
emphasize here that for this quenching scheme [25|, the 
exponent zi is fixed (=3) irrespective of the value of a. 

Since in the limit of r — > 00, the maximum contri- 
bution to the defect comes from the vicinity of ^q, we 
linearize the Hamiltonian (6) [l!] around t ^ to 



t_to. 

Toff 



(7) 



with Toff = Tfc-2("-i)/"/a and r^ff = r for a = 1. The 
higher order terms lead to a faster decay of the defect 
density and hence can be dropped in the limit of large r. 
The excitation probability is now readily obtained using 
the Landau-Zcner formula pk = exp(— Trfc^Tcg) which is 
integrated over all modes k to derive the scaling of the 
defect density given as 



dfcpfc ^ (I)-"/[6-+2(l-a)] ^ (I)-o/[2(2a+l)]^ (g) 



It is straight forward to generalize to a d-dimensional 
quantum MCP with viZi ~ 1 described the Hamiltonian 



(l-rsgn(t) 

T 



k'')&'' + |t/T|"sgn(t)fc^c 



(9) 
(10) 



where the MCP is reached at time t = due to the si- 
multaneous variation of two parameters Ai and A2 under 
the quenching scheme Ai = A2 = |i/r|"sgn(t). Compar- 
ison with the Hamiltonian (5) gives Ai=/i— 1, A2=7, 
zi ~ 2 and /3 = 1. Linearizing around to and using 
Toff = r/s^^i("^-'^)/"/a, we similarly get 



„ ^ ^-ad/[2z2a+zi(l-a)]. 



Although scaling given in Eq. (11) matches identically 
with the scaling relation Eq. (4) of reference [l^ for a = 1 
(with do = and viZi = V2Z2 ~ 1), there is a small 
difference for a ^ 1. This mismatch stems from the 
fact that we have used linearization of the Schrodingcr 
equations in the vicinity of to (or Tcff) which necessarily 
yields a term zi(l — a) in the KZS for a 7^ 1. We do 
also note that for a < 1, Toff — )• in the limit — > 
rendering the process of linearization inappropriate. In 
the following, we shall restrict our analytical calculations 
to a > 1. 
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FIG. 3: (color online) The defect density obtained by numer- 
ical integration of the Schrodinger equation for path a = 1.5 
(slope ~ 0.29, green line), a — 2 (slope ~ 0.38, red line) and 
a = 2.5 (slope ~ 0.39, blue line). In the inset, we show a 
comparison between numerically obtained (dashed line) and 
analytical obtained (solid line) slopes for different a. The 
numerical results are in good agreement with the analytical 
predictions; the observed small mismatch is due to neglecting 
the higher order corrections in the process of linearizing. 
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FIG. 4: Defect density obtained by numerical integration of 
the Schrodinger equation for special path III; a = 1 (slope 
^ 0.73), a = 1.2 (slope « 0.63) and a = 1.5 (slope 0.5). 
Inset shows the comparison between numerically obtained and 
analytical obtained slopes for different a. The numerical re- 
sults are in good agreement with the analytical predictions. 



Now for our quenching scheme (2), the Hamiltonian 
(3) takes the form 



(11) Hk = (l-rsgn(t) + k')a'- + -ka\ 



(12) 



We note that the energy gap of Hamiltonian (12) is 
minimum at a quasicritical point for t = to where 
|to/T|"sgn(t) fc2 = only if a < 2 while for a > 2, the 
minimum gap occurs right at the MCP (i.e., at t = 0) 
(see figures (1) and (2)). Linearizing around t = to, we 
get 



Hk = 



t-tp 

Toff 



(13) 



with Teflf = rfc ^^1°" ja. Here, we have ignored the 
time dependence of the off-diagonal term in favor of the 
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minimum value at t = to scaling as k'^/°'~^^ which is 
appropriate for 1 < a < 2 but not for a > 2 since 
Teff/(r/fc) - A:(2-")/" ^ as /c ^ if a < 2. This 
implies that Toff dictates the scaling behavior of the de- 
fect for the dynamics of Eq. (12) for 1 < a < 2 while for 
a > 2 a crossover to the scaling behavior dictated by kt/r 
is expected. Further, the exponent Z2 = ^2(0) = 2/a+ 1 
depends on a though zi is completely independent of the 
path and fixed (=2). We emphasize that the existence 
of the path dependent Z2 is the crux of the quenching 
scheme (2) which leads to striking consequences shown 
below. We already note that a = 2 is a marginal case 
where Z2 becomes equal to zi 

Application of the Landau-Zener formula to Eq. (13) 
similarly yields a scaling relation for the defect density 
(71,2) for the quenching scheme (2) given by 



712 = / dfcexp(— 7rfc°+ Toff) ~ ( — ) 

a 



(14) 



so that the exponent varies continuously with the pa- 
rameter a. In the limit a = 1, we recover the well known 
-1/6 in while for a > 2, we get " ^ ^-i/3 



scaling n < 

p^ . [14 1 the system hits the MCP vertically and does no 
longer cross any quasicritical point along the path and 
hence the quenching scheme boils down to driving the 
system across the MCP along the Ising critical line (Tsj 
when the term It/rj" in Eq. (12) can essentially be ig- 
nored. The defect density as a function t for different a 
obtained by numerical diagonalization of the Schrodinger 
equations in case of full path is shown in Fig. 3, and 
the results are in good agreement with the prediction of 
Eq. (14). We note that the linearization process retains 
only the slowest r-dependence of ti, while contributions 
from the higher order terms decay must faster. However, 
this can be numerically observed only in the limit of suf- 
ficiently large r. We reiterate that the scaling relation 
given in (14) is inappropriate for a < 1 though our nu- 
merical results do indicate the existence of continuously 
varying exponent that decreases with decreasing a. 

Generalizing to the d-dimensional Hamiltonian (9) 
with Toff = rfc~^i'^"~-'^^/"/a, we get the relation 



"2 



-aci/[2z2a+zi(l-a)] 



(15) 



Although relation (15) closely resembles (11), it is to be 
noted that in our quenching scheme the exponent 22 (= 
zi/a -I- /?) and hence, the scaling exponent depends on a. 
up to a = Qfc when 2;i(ac) = ziipcf^ i.e., ac = zi/(zi — /3). 
Finally, we consider the quenching scheme 



h{t) = 1 + I7I" and 7 = t/r. 



(16) 



in which the system lies in the paramagnetic phase 
throughout the quenching process and touches the MCP 
at t = (see path III in Fig.l). The case a = 1 was 
studied in reference [2^ and a scaling relation of the 
form n ^ r~^/^ was obtained. This anomalous scaling 
behavior was justified using the argument of a dynam- 
ical shift of the center of the impulse region which 



makes it asymmetric as opposed to the conventional KZS 
such that non-critical energy states contribute to the de- 
fect with pk ^ fc''o(") where do{a) is the additional effec- 
tive dimension which eventually appears in the KZS. In a 
similar spirit, we use the the LZ formula for Hamil- 
tonian (13) for the half path so that the singularity at 
t = is avoided; this provides the correct scaling for 
the full path since the symmetry of the initial and final 
states are the same. Expanding the transition proba- 
bility (Eq. (7) of HI) in Taylor series around Tj = 0, 
we get dominant contribution given as pk ~ I'^/T/P 
where the dimensionless time T/ = (}? ^ t) j ^ t^s and 
the dimensionless coupling = ^(2/0+1)^^^^^ Using 
Teff = Tfc-2("-i)/"/Q,, we get pfc = (l/a2)fc2(2-")/". The 
approximation is valid up to fc = fcmax for which Tf = 1 
given by fcmax ~ r""/^*-""*"^-*. The scaling of the defect 
density (rigp) for these special paths is hence given by 

risp = / dkpk ^ fcmSx ~ T 2<i+°) , (17) 

where in the limit a = 1, we recover the scaling t~^/^ as 
in [2^. We conclude that for quenching scheme (16) the 
dimensional shift do^a) = 2(2 — a) /a, varies continuously 
with a; do{a) = 2 for a = 1 and vanishes for a > ac{= 2) 
when once again the impulse region becomes symmetri- 
cal. Interestingly, for these special paths also we obtain 
a continuously varying exponent which saturates to 1/3 
in the limit a — > ac. Fig. 4 shows the defect density as 
a function of r for different a for the special path HI, as 
obtained from numerical integration. 

We shall now derive the generalized scaling relation 
for a quantum system with viZi ^ 1 using adiabatic per- 
turbation theory We note from Eq. (13) that the dy- 
namics across a MCP can be described as a linear quench 
with a different Tgff for each k mode across a quasicritical 
point which provides the most dominant contribution to 
the defect. 

For a general Hamiltonian parametrized by A = t/r, 
we can write the wave function as ■(/; = X^n '^p(0'/'('^)- 
the above basis the Schrodinger equation oecomes 



i dap{t) i 
T dX T 



(18) 



where Ep{X) is the eigenenergy of H{X). A unitary trans- 
formation 



ap(t)=ap(A)e-*^/'^-(^)^^ 
transforms eq. (18) to 
ddp{X) 



(19) 



dX 



= - y a„ ( A) (p| |<7) e'^- .(20) 
^-^ dX 



Now taking into account that initially the system was in 
the ground state, which makes a single term dominate 
the sum in eq. (20), and expressing the above equations 
in the momentum basis, we get for our case. 



d'^k 



{2n) 



dX{k\±\0)e^-^''^\^\-'-)r'ix'SE,ix') 
oX 



.(21) 



5 



where near a quasicritical point, SEj:{\) ~ 

|A|^^'^^F'(|fc|^VI^I''"'')- Here A denotes the devi- 
ation from the quasicritical point and F'{x) ^ x 
for large x. Again, we have (fc|^|0) jGik/X"^) 
where G(0) is a constant. Using these relations, 
defining A' = X/\k\'^^'^^ and eventually rescaling 
\k\ \k\T°"^2/Mz2^2+i}+zi,^2{i-a)] ^ get 

n2 ^ ^-dai'2/[a(22i'2 + l)+2l!^2(l-a)] (^22) 

which reduces to Eq. (f5) for 1^222 = 1- Of course, the 
scaling relation (22) holds up to a value of a = such 
that zi{ac) = 22(0^0) and for a > ac gets modified to 
n ~ ^-da^ui/la^+uizi]^ For special paths (e.g., path III), 
we note that the minimum gap is given by exponents as- 
sociated with the MCP and derive the generalized scaling 
relation given by 

jl^^ ^ ^-(li+do(Q))Ql'l/[Q2ll/i + l] ^2^"j 

where we have included the dimensional shift do (a) which 
varies continuously with a for a < ac, while for a > ac 
the scaling saturates to n T--aaiyid/{a,^izi+i) ^ 



The main results of the paper are summarized in 
Eqs. (14,15,17) and (22-23). In conclusion, we propose 
a different quenching scheme across the MCP and de- 
rive a KZS with an exponent continuously varying with 
the path as a consequence of path dependent effective 
dynamical exponents associated with the quasicritical 
points. This continuously varying exponent up to a crit- 
ical value of a is the cardinal feature of this paper. Also, 
our studies establish that unlike the conventional KZS, 
the path-dependent effective exponents which usually do 
not appear in studies of the equlibrium quantum critical 
behavior, do appear in the scaling relation of the defect 
density in special situations like quenching across a mul- 
ticritical point. 
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